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Abstract 

A theory is proposed to describe the competition among antiferromag- 
netism (AF), spin glass (SG) and Kondo effect. The model describes two 
Kondo sublattices with an intrasite Kondo interaction strength Jk and 
an interlattice quantum Ising interaction in the presence of a transverse 
field r. The interlattice coupling is a random Gaussian distributed vari- 
able (with average —2Jo/N and variance 32J^/N) while the T field is 
introduced as a quantum mechanism to produce spin flipping. The path 
integral formalism is used to study this fermionic problem where the spin 
operators are represented by bilinear combinations of Grassmann fields. 
The disorder is treated within the framework of the replica trick. The 
free energy and the order parameters of the problem are obtained by us- 
ing the static ansatz and by choosing both Jo/ J and F/J ~ (Jk/J)^ to 
allow, as previously, a better comparison with the experimental findings. 
The results indicate the presence of a SG solution at low Jk/J and for 
temperature T < Tf (Tf is the freezing temperature). When Jk/J is 
increased, a mixed phase AF-I-SG appears, then an AF solution and fi- 
nally a Kondo state is obtained for high values of Jk/J. Moreover, the 
behaviors of the freezing and Neel temperatures are also affected by the 
relationship between Jk and the transverse field F. The first one presents 
a slight decrease while the second one decreases towards a Quantum Crit- 
ical Point (QCP). The obtained phase diagram has the same sequence as 
the experimental one for Ce2Au\-xCoxSi-i, if Jk is assumed to increase 
with X, and in addition, it also shows a qualitative agreement concerning 
the behavior of the freezing and the Neel temperatures. 

1 Introduction 

The competition between RKKY interaction and Kondo effect is recognized as 
the fundamental mechanism to explain several properties of Ce and U com- 
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pounds. Nevertheless, there are also strong experimental evidences that the 
presence of disorder in those compounds can be source of quite novel effects 
mil 121 EE]- For instance, deviations of the Fermi-Liquid (FL) behavior are 
observed in the magnetic, thermodynamic and transport properties. Recently, 
theories have been proposed to explain how disorder can produce such devi- 
ations [HIE]. However, it is still not clear whether disorder itself can be the 
origin of such deviations. Other prospective has adopted the proximity to a 
QCP jS] as a possible origin of the Non-Fermi Liquid (NFL) behavior. Besides 
this complicate scenario with NFL behavior, some of these physical systems 
present a highly non-trivial manifestation of disorder which is frustration. As 
consequence, there are many phase diagrams with ferromagnetism (FM), an- 
tiferromagnetism (AF), spin glass (SG) and a Kondo state with the complete 
screening of the localized spins. This raises the question of whether it is possible 
to build up theories able to account for those phase diagrams. 

The experimental situation of Ce and U alloys which present together a spin 
glass phase, an antiferromagnetic one and finally a Kondo one is complicate, 
especially when the phase changes occur with changing the matrix composition 
of the alloy. Moreover, there is certainly an important local effect in the spin 
glass, as recently evidenced by local measurements in CeCui^xNix alloys [5]. In 
particular, there is a change from AF to Kondo and then to SG with increasing 
X in UCu^-xPdx jlOj and a different situation in Ce2Aui^xCoxSi3 alloys 
More precisely, in the case of Ce2Aui-xCoxSiz alloys, the glassy behaviour 
is not favored when the chemical induced disorder is initially increased. In the 
experimental phase diagram, for < a; < 0.45 and low temperature, a spin glass- 
like state appears. When the Co doping is increasead, the system experiments a 
transition to an AF phase. In the interval 0.45 < a; < 0.9, the Neel temperature 
Tjv decreases until reaching a QCP at a; = 0.9 with no trace of NFL behavior. 
For X > 0.9, the magnetic moments are supressed due to the Kondo effect. 

Recently, a theoretical approach has attempted to describe the interplay 
between AF and SG in a Kondo lattice using the same framework previously 
introduced to study the SG in the Kondo lattice ^^j- The proposed model has 
been a two-Kondo sublattice with an intrasite exchange interation with strength 
Jk and a random Gaussian Jij intersite interaction only between distinct sub- 
lattices. Two quite important points are introduced in this approach. The first 
one is that there is no hopping of conduction electrons between distinct sub- 
lattices. Therefore, the antiferromagnetic solutions in this model are entirely 
associated to the coupling among localized magnetic moments. The second one 
takes the relationship Jq/Joc (Jk/J)'^, because the two exchange terms cannot 
be considered as completely independent from each other 16 . As a consequence, 
when the strength of the intrasite exchange interaction Jk increases (in units 
of J), the degree of frustration J/Jq decreases. This last parameter controls 
the emergence of SG or AF solutions in the problem. The results have been 
shown in a phase diagram T/J (T is the temperature) versus Jk/ J- For small 
Jk/J and low temperature, there is a SG phase. When Jk/J is increased, an 
AF solution has been found. Finally, for Jk/J > 15, the Kondo state becomes 
dominant. Thus, we can conclude from the previous calculation |12| that the 
sequence of phases SG, AF and Kondo mimics qualitatively the experimental 
phase diagram of Ce2Aui-xCoxSi3 if Jk/J is associated with the content of 
Co. 

But, we can say that the previous agreement was only qualitative to provide 



2 



the correct sequence of phases for the alloys Ce2Aui-xCoxSi3 JT], because 
in the previous approach there is no mechanism able to produce a QCP. As 
a consequence, this approach cannot provide the proper behavior of the Neel 
temperature which should decrease towards a QCP. Thus, the purpose of the 
present paper is to obtain specifically a quantitative description of the phase 
diagram of Ce2Aui-xCoxSi3 alloys and in particular of the QCP, by adding a 
quantum mechanism to produce spin flipping, given by an additional transverse 
field r in|. 

We take here the replica symmetry mean-field approximation and we use the 
static approximation. Indeed, such an approximation is subject to criticism, but 
it was shown in many previous works that it yields a good description of the 
competition between the Kondo effect and magnetism in the Kondo lattice mod- 
els. Moreover, the quantum Ising spin glass in a transverse field has the same 
critical behavior as the M-component quantum rotor which is exactly solvable 
in the M — oo limit JS]- In this model, the line transition is given by the 
singularity of the zero-frequency that is equivalent to the static approximation. 
The static approximation is, therefore, relatively well justified for the study of 
phase boundaries, which is really the purpose of our paper in the specific case 
of the SG-AF-Kondo competition in alloys such as Ce2Aui^xCoxSi3. 

Then, we will assume a relationship T cx J^; this choice, which has been 
already used in references ^\ E], is taken to account for the fact that the 
intrasite exchange interaction is able to produce both the Kondo effect and the 
RKKY interaction. The transverse field plays a role similar to the spin flipping 
part of the Heisenberg model. One of the main achievements of this approach 
has been to show the presence of a QCP in the competition between SG and 
the Kondo state when P is enhanced and to improve, therefore, considerably 
the description of experimental results. 

In this paper, we study the SG-AF competition in the disordered Kondo 
lattice, described by the hopping of the conduction electrons only inside each 
sublattice, an intersite exchange interaction, given by an Ising-like term only 
between different sublattices and a transverse field P applied in the x-direction, 
in order to obtain a mechanism able to produce a QCP. We use also the static 
approximation and the relationships Jo cx and P oc already introduced in 
references j^lEl- Therefore, the initial many parameters problem is reduced 
to a one parameter which allows a better comparison with experimental results. 

The comparison of our theoretical results to experimental data is delicate 
for two reasons. First, as we have just discussed, we take a J^-dependence of 
both Jo and P, which is clearly a crude approximation but which has already 
given interesting physical results in reference p~61| . Then, the second question 
concerns the fact that Jk is assumed to increase with the concentration of cobalt 
in Ce2Aui-xCoxSi3 alloys, as previously done with the concentration of Nickel 
in CeCui-xNix alloys. Indeed, the role of the so-called "chemical pressure" is 
less clear than that of the regular pressure, but obviously the spin glass phase 
originates in such alloys only from the disorder of the matrix and we are obliged 
to make an assumption of the relationship between Jk and the relative host 
concentration x, which was previously successful in the case of CeCui^xNix 
alloys reference il7. . This question is in fact not easy to answer theoretically. 

In section 2, we introduce the model and the relevant order parameters for 
the problem. The saddle point free energy in terms of the order parameters is 
obtained so that it allows to obtain a phase diagram giving temperature T/J 
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versus Jk/J- The results and the final discussion are presented in the last 
section. 



2 General Formulation 

We have considered in this work a model given by two Kondo sublattices A and 
B with a random coupling J^- only between localized spins in distinct sublattices 
|12[I19| . There is also a transverse field F coupled with the localized spins in both 
sublattices As mentioned in section 1, the hopping of conduction electrons 
between two different sublattices is not allowed for simplicity. Therefore, the 
Hamiltonian is: 



p=A,B i,j (T=ti i 



(1) 



where i and j sums are over A'' sites of each sublattice. The intersite coupling 
Jij is a random variable following a gaussian distribution with average —2Jq/N 
and variance 32J'^/N. The spin operators present in Eq. are defined as 



Ref. [n]: 



'i,p 



fi,p,]'fi:P^l h,p,i fi,P,i 



and S'f „ = i 



fi,p,^fi.Pd ~^fi,p,l fi,PA 



^i.p,a 



where di p 

ifi p ai fi,p-<^) are the creation and destruction operators for conduction 
(localized) fermions. 

The partition function is given in the Lagrangian path integral formalism 
in terms of Grassman variables |12l I13| '0i,p,o-(T) for the localized fermions and 
'4^i,p,ij{T) for the conduction ones. The Kondo interaction is treated following 
Ref. 112]. The Kondo state is described by complex order parameters Ap.o- — 
'h'^i (T ^ fi p cr''Pi,p,cr > which arc introduced in the partition function using 
the integral representation of delta function. Therefore, the partition function 
when Xnrr « A„ becomes: 



= cxp{-2Npjki\XA\^ + I Asp)} f Yl Di^p^p) exp 

A D 



A{stat) 



P=A,B 



(2) 



where Zq is the partition fmiction of the free conduction electrons. 



5..(w) 



e,-(^) + /3^j.,5f,^^,s, 



ej(w) = [ V'r,yl,T(^) ^lAd'y^) ^Iba{^) i^lB, li^) 



(3) 
(4) 



and 



is a 4 X 4 matrix given by 





-1 











=B,ij ^ ' 



(5) 
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with 

(6) 

The notations / and represent the identity and the null 2x2 matrices, re- 
spectively. 

The free energy is given by the replica method. The procedures are quite 
close to reference |12) and in particular to its appendix, but with now the pres- 



ence of the transverse field T which modifies the matrix 



Therefore, 



the averaged free energy within the replica symmetric theory can be obtained 
as 



PF = I3Jk{\Xa\'' + IAbH - P'^ilAqB - qAqs) - P^ruAmB - lim^ 



1 



N 



X In / n Di.,ADU,B n / Dzl^Dzls exp 



n^o 2Nn 
1 



^Indet gAuj\hl^) 



(7) 



where the matrix 



G..(w|ft-f ) in Eq. Q is the same matrix given in Eq. 
(O, except that the matrix F. . (cj) (see Eq. (jSjl) is replaced by F.. {lo) + 

a/3J/i^/yI(a-+(T),-a)). InEq. ^ Dy ^ dy e^-vV^) / iv ^ L z) ■ 
The internal field /if ^ applied in the sublattice p ^ A, B \s defined as 



h°' 

2,P 



(8) 



In Eq. ||SJ), Qpi, Xp' = PXp' ^^'^ ™p' &re the SG order parameter, the static 
susceptibility and the magnetization, respectively. It should be noticed that the 
internal field applied in sublattice p depends on the behavior of qp>, Xp' and 
uipi with p' 7^ p. That is a direct consequence of the choice made in the model 
where only interlattice frustration has been considered in the present work. 

r 1 ~^ 

Finally, to solve Eq. |7J, the matrix G_. .(uj\hf ) is substituted by 

G^^{lu\ hfp)] ^. In fact, this approximation |12l I13L ITT) represents a Kondo 

sublattice p in the presence of a constant magnetic field /if^ where the effects 
of the Kondo sublattice p' are placed. Therefore, in Eq. (0, we have: 



Indet 



a,,HKp)Y = NT.^'^'^'^ia.MKp) 



(9) 



The sum over the Matsubara's frequencies in Eq. (0 can be now performed 
as usual ^1 El E| ■ The proposed decoupling also allows us to use Fourier 
transformation which can be evaluated by assuming a constant density of states 
for the conduction electron band p{e) = 1/2_D for —D < e < D. Thereby, the 
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free energy is obtained as: 



PJo 1 
—rriAmB - ^ 



D^3s In 



n 

p=A,B 



where 



with 



E{Hp) 



1 

Jd j_ 



dxln{cosh( P ) + cosh(VA)} 

fiD 2 



1 



(10) 



(11) 



(12) 



and -ffp = + /i^^ with the internal field /ip given by Eq. © (without here 

the site and replica indices). The saddle point order parameters follow directly 
from equations (|1()|I -H12 |I . 

The limit of stability for the order parameters with replica symmetry is given 
when the Almeida- Thouless eigenvalue \at becomes negative |14l I15| : 



/oo 



V=A,B ' 



-\ 2 



(13) 



where 



J— oo J —oo orip orip J ~oo 



E(H„) dE{Hp) 2 
dhp ^ ■ 
(14) 



3 Results and Conclusions 

The central issue of the present work is to account for the general features of the 
experimental phase diagram of Ce alloys such as Ce2Aui-xCoxSi3 11 which 
present with increasing Co concentration the successive sequence of the SG 
phase and then the AF phase with the Neel temperature Tjv decreasing towards 
a QCP. In principle, the solutions for the saddle point order parameters qp, Xp, 
rup and Ap (p = A, B) should be found in a parameter space with axes given 
by Jo/ J, Jk/J and T/J (J and D are kept constant). 

Therefore, in order to identify the role of each parameter in the problem, 
the solutions for qp, Xp, rUp and Ap {p = A, B) are shown in figure 1 as a phase 
diagram T/J versus Jk/J for several values of Jq/ J and T/J. The SG solution 
corresponds to qa = qs 7^ (with tua = tob — 0) while the AF solution is 
equivalent of ttia = —mB ^ and qA — qs Y 0- The Kondo state in the 
phase diagram corresponds to the situation where the two Ap {p = A, B) are 
different from zero. In figure l.a, Jo = and F = or F/J = 1.5, this situation 
corresponds to the studied case in Ref. in which there is no AF solution. 
For F = 0, the phase diagram displays three solutions. At high temperature and 
small Jk / J , a paramagnetic (PARA) solution is found. When the temperature 
is decreased, there is a transition to the SG phase at Tg, which coincides with 
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Figure 1: Theoretical phase diagram T/J versus Jk/J for several values of 
Jo/ J: (a) Jo/ J = 0, (b) Jo/ J = 1.5, (c) Jo/ J = 1.8 and (d) Jq/J = 2.0. The 
full lines correspond to the results obtained for T/J = 0.0, while the dashed 
ones are for T/J= 1.5. 
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the Tf temperature that locates the Almeida-Thouless hne (the temperature 
where the Almeida-Thouless eigenvalue is zero). If Jk/J is enhanced, there 
is a transition at Jk = Jki{T/J) to a region where the Kondo solution is 
dominant. When F/J is tuned on, there is a clear effect on the transition lines 
Tg and Jki{T/J), the first one is decreased while the second is displaced to a 
larger Jk/J value. 

When Jo/ J is enhanced, as it is shown in figures l.b - l.d, an AF solution for 
the order parameters appears below the Neel temperature T/v for small Jk/J- 
For this particular situation, the set of temperatures Tjv, Tf and Tg is function 
of F/J and Jq/J. For F = 0, the temperature Tjv increases with Jq/J while 
the replica symmetric transition temperature Tg and the Almeida-Thouless line 
Tf decrease with Jq/J. Nevertheless, similar to Ref. [131, ^/ becomes larger 
than Tg. Therefore, in such circumstances, there is the onset of a mixed phase 
AF-I-SG at Tf. This mixed phase represents a replica symmetry breaking SG 
phase with non-zero spontaneous sublattice magnetization. Finally, at Tg there 
is a transition to a pure SG phase. For increasing Jk/J, a Kondo solution at 
Jk > Jki{T/J) is found again. When the transverse field F/J is tuned on, its 
effect is basically to decrease simultaneously Tjv, Tf and Tg and displace the 
transition line Jki{T/J) as before. 

However, to compare the results obtained in this work with the experimental 
results, it is better to reduce the number of parameters which are in fact not 
independent from each other, because both the Kondo effect and the RKKY 
can be generated directly from the intrasite exchange interaction. Thus, we 
choose a relationship among the parameters Jo/ J, Jk/J and F/J, similar to 
that used in Ref. |E1 where a relationship Jq/J oc [Jk/JY has successfully 
reproduced the experimental sequence of phases. Nevertheless, in Ref. |12| . 
we cannot explain at all the decrease of the Neel temperature T^r towards a 
QCP. In the present work, the transverse field F has been introduced to provide 
a simple mechanism able to reproduce the spin flipping part of the Heisenberg 
model (ITj. Therefore, it would be also natural to assume F/J oc {Jk/ JY in the 
present problem. It should be remarked that this choice introduces an additional 
complexity. It has been shown that F/J and Jo / J enforce different effects in a 
previously studied competition between AF and SG. Within a certain range of 
Jo/ J, this new choice can favor SG, SG+AF or AF, while F/J tends to destroy 
the three phases leading the transition temperatures to zero |20| . 

The solution of the set gp, Xj "^p and \p with Jo/J oc F/J c>c {Jk/JY is 
shown in figure 2.b giving the phase diagram T / J versus Jk / J ■ In particular, 
the Neel temperature Tjv = T^^^Jk) can be computed expanding the saddle 
point equations in powers of nip. At the second order critical line Tat, we can 
make Qp = Xp ~ 0, therefore X = Xp = ^/{PcJq) and 



where (3c = 1/JW and m = niA = —ms. For Tjy tending to zero, the critical 
value of Fc is given by analytical solution of Eq. H15(l : Fc = Jq + 2J'^ / Jq. For 
< Jk/J < 13.5, figure 2.b shows the presence of the SG solution in which Tg 
and Tf coincide. Then, in a very small interval of Jk/J, Tf ^ Tg. Therefore, 
the mixed phase AF-I-SG appears after the SG phase. Finally, in the interval 
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Figure 2: Comparison between experimental and theoretical phase diagrams: 
the upper figure is an experimental phase diagram of (762^^1-0:^^02; 5*1 3 alloys 
|11| : the lower figure is a theoretical phase diagram T/J versus Jk/J for the 
relations V = 0.011( J^/ J)^ and Jq/J = 0.73r/J, where the dotted lines are 
the extrapolations carried for lower temperature. 
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13.9 < Jk/J < 17.0, only the AF solution is present. When Jk/J increases, 
the Neel temperature decreases until a QCP. Above this point onwards, there 
is a transition line to the Kondo state characterized by = As ^ JTj. At 
high temperature, the solutions of the order parameters are qp = 0, rup = and 

Xp = Oip^A, B). 

The situation described in the present replica symmetry mean field theory 
can be analysed as follows: for a low Jk/J value (corresponding to high value 
of J/ Jo), frustration is dominant and the long range internal field is given by 
the SG component of hp (see Eq. ©). When Jk/J is increased, the degree of 
frustration J/ Jq is decreased and, therefore, the localized spin of the sublattice 
A (B) starts to couple with an internal field which depends on the negative 
magnetization of the sublattice B {A) in the AF ordering y[2i- However, this 
process also implies the increase of the transverse field F and consequently the 
freezing temperature (Tg) is slightly decreasing. On the other hand, the Neel 
temperature (Tjv) is deeply affected and decreases towards a QCP. Finally, for 
enough high Jk/J, there is complete dominance of the Kondo effect ^21E1E|- 
Finally, as we have already explained in the introduction, the comparison be- 
tween theory and experiment is more delicated in the case of disordered alloys 
where the different phases occur as a function of the relative matrix concen- 
tration than in the classical case where increasing pressure makes Jk increase. 
But here we consider that Jk increases with increasing Co concentration , as 
we have successfully assumed that Jk increases with increasing Ni concentra- 
tion in the case of CeCui-xNix alloys reference |17| . Under this assumption, 
the theoretical results given in figure 2.b agree quite well with the experimental 
phase diagram of Ce2Aui^xCoxSi-i alloys (see figure 2. a) and in particular 
we obtain the same order for the sequence of phases and, mainly, the correct 
behavior of the Neel temperature. 

To conclude, we have used a two-sublattice model with two exchange inter- 
actions, an intrasite exchange and a random intersite between localized Ising 
spins in the presence of transverse field. The main goal of the present work has 
been to reproduce some fundamental aspects of the phase boundary contained 
in the experimental phase diagram of the Ce2Aui^xCoxSi^ alloys jll| . Finally, 
the theoretical results shown in figure 2.b account quite well for the most im- 
portant part of the experimental phase diagram of Ce2Aui-xCoxSi'i alloys |11| 
with the cobalt concentration (figure 2. a) and the agreement between the two 
figures shows a clear improvement performed by the present model with respect 
to previous ones. In fact, important questions are not really solved, such as 
for example the different possible Kondo-AF-SG sequences found in different 
disordered alloys, the relationship between the parameters of the model and the 
matrix concentration which is certainly less clear than the pressure dependence 
and finally the precise local nature of the spin glass phase. Further work will 
be, therefore, necessary to improve the theoretical description and to obtain 
new phase diagrams of Ce or C/ alloys or compounds with either the matrix 
concentration or the external pressure. 
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